AFRL-OSR-VA-IR- 2015-0220 


Nonlinear  dynamics  of  self- pulsing  all-solid-state  lasers 


Alejandro  Hnilo 

UBAIEC  SA. 


07/ 06!/ 2015 
Final  Repoit 


DISIRIBUnON  A:  Distribution  approved  for  public  lelease. 


AirFbrce  Research  Laboratory 
AF Office  Of  Scientific  Research  {AF03R)/  lOS 
Ariington,  Viiginia  22203 
AirFbrce  Materiel  Command 


Report  Documentation  Page 


Form  Approved 
0MB  No.  0704-0188 


Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  0MB  control  number. 


1.  REPORT  DATE 

06  JUL  2015 


2.  REPORT  TYPE 


4.  TITLE  AND  SUBTITLE 

Subtitle:  Nonlinear  dynamics  of  self-pulsing  all-solid-state  lasers 


6.  AUTHOR(S) 


7.  PEREORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Centro  de  Investigaciones  en  Laseres,y  Aplicaciones  (CEILAP), , , 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 


3.  DATES  COVERED 

27-03-2013  to  31-03-2015 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


10.  SPONSOR/MONITOR’S  ACRONYM(S) 

11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

The  Project  involved  the  study,  both  experimental  and  theoretical,  of  extreme  events  (EE,  sometimes 
named  Optical  RogueWaves)  in  two  types  of  self-pulsing  lasers:  1)  all-solid-state  (i.e.,  diode-pumped) 
Nd:YV04-FCr:YAG  (slow  saturable  absorber);  and  2)  Kerr-lens-mode  locked  Ti  Sapphire  (fast  saturable 
absorber).  The  main  tool  is  the  analysis  of  the  time  series  of  the  laser  pulses  and  of  the  images  of  the  laser 
spots.  Regarding  (1),  observations  on  a  specially  designed  and  constructed  prototype  showed  the  existence 
of  EE  in  the  pulse  energy  and  also  in  the  time  separation  between  pulses.  The  two  types  of  EE  are 
uncorrelated:  a  high  energy  pulse  is  not  preceded  by  a  longer  pumping  time.  Therefore,  the  problem  is  to 
identify  the  ???reservoir???  for  the  energy-EE.  Correlations  were  identified  between  the  spatial 
complexity  of  the  laser  spot,  the  Eresnel  number  of  the  cavity,  the  dimension  of  embedding  of  the  time 
series  and  the  appearance  of  EE. 

15.  SUBJECT  TERMS 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 

18.  NUMBER 

19a.  NAME  OE 

ABSTRACT 

OF  PAGES 

RESPONSIBLE  PERSON 

a.  REPORT 

unclassified 

b.  ABSTRACT 

unclassified 

c.  THIS  PAGE 

unclassified 

Same  as 
Report  (SAR) 

20 

Standard  Form  298  (Rev.  8-98} 

Prescribed  by  ANSI  Std  Z39-18 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
0MB  No.  0704-0188 


The  public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including 
suggestions  for  reducing  the  burden,  to  the  Department  of  Defense,  Executive  Service  Directorate  (0704-0188).  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no 
person  shall  be  subject  to  any  penalty  for  falling  to  comply  with  a  collection  of  Information  If  It  does  not  display  a  currently  valid  OMB  control  number. 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ORGANIZATION. 


2.  REPORT  TYPE 

Final  Technical  Report 


1.  REPORT  DATE  ('DD-MM-YV'yy; 

12-05-2015 


4.  TITLE  AND  SUBTITLE 

Title:  Final  Report  2015 

Subtitle:  Nonlinear  dynamics  of  self-pulsing  all-solid-state  lasers 


3.  DATES  COVERED  (From  -  To) 

27  March  2013-31  March  2015 


5a.  CONTRACT  NUMBER 


5b.  GRANT  NUMBER 

FA9550-13-1-0120 


5c.  PROGRAM  ELEMENT  NUMBER 


6.  AUTHOR(S) 

Alejandro  A.  FInilo  and  Marcelo  G.  Kovalsky. 


5d.  PROJECT  NUMBER 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

1 .  Centro  de  Investigaciones  en  Laseres  y  Aplicaciones  (CEILAP). 

J.B.  de  La  Salle  4397,  (B1603ALO)  Villa  Martelli,  Argentina 

2.  UBATEC 

Av.Roque  Saenz  Pena  938  6to  piso,  (Cl 035 AAR)  Buenos  Aires,  Argentina 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Air  Force  Office  of  Scientific  Reseach 

875  North  Randolph  St.,  Suite  325,  Room  3112;  Arlington,  VA  22203-1768,  USA. 


8.  PERFORMING  ORGANIZATION 
REPORT NUMBER 


1 0.  SPONSOR/MONITOR'S  ACRON  YM(S) 

AFRL/AFOSR 


1 1 .  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 


12.  DISTRIBUTION/ AVAILABILITY  STATEMENT 

DISTRIBUTION  A  APPROVED  FOR  PUBLIC  RELEASE 


14.  ABSTRACT 

The  Project  involved  the  study,  both  experimental  and  theoretical,  of  extreme  events  (EE,  sometimes  named  Optical  Rogue  Waves)  in  two  types 
of  self-pulsing  lasers:  1)  all-solid-state  (i.e.,  diode-pumped)  Nd:  YV04+Cr:YAG  (“slow”  saturable  absorber);  and  2)  Kerr-lens-mode  locked  Ti: 
Sapphire  (“fasf  ’  saturable  absorber).  The  main  tool  is  the  analysis  of  the  time  series  of  the  laser  pulses  and  of  the  images  of  the  laser  spots. 

Regarding  (1),  observations  on  a  specially  designed  and  constructed  prototype  showed  the  existence  of  EE  in  the  pulse  energy  and  also  in  the  time 
separation  between  pulses.  The  two  types  of  EE  are  uncorrelated:  a  high  energy  pulse  is  not  preceded  by  a  longer  pumping  time.  Therefore,  the 
problem  is  to  identify  the  “reservoir”  for  the  energy-EE.  Correlations  were  identified  between  the  spatial  complexity  of  the  laser  spot,  the  Fresnel 
number  of  the  cavity,  the  dimension  of  embedding  of  the  time  series  and  the  appearance  of  EE. 

Regarding  (2),  two  coexistent  modes  of  operation  are  observed:  transform- limited  pulses  (PI)  and  chirped  pulses  (P2).  EE  are  observed  only  in 
P2.  We  found  that  EE  arise  after  a  threshold  similar  to  the  modulational  instability  is  crossed.  EE  exist  in  P2  only,  because  PI  is  unstable  for  the  par 


15.  SUBJECT  TERMS 

Final  Technical  Report  for  the  period  2013-2015  of  the  project  "Nonlinear  dynamics  of  self-pulsing  all-solid-state  lasers".  Study  of  optical  rogue 
waves  in  two  types  of  self  pulsing  lasers:  diode-pumped  Nd:YV04+Cr:YAG  and  Kerr-lens  mode-locked  Ti:Sapphire. 


16.  SECURITY  CLASSIFICATION  OF: 

a.  REPORT 

b.  ABSTRACT 

c.  THIS  PAGE 

u 

u 

u 

17.  LIMITATION  OF 
ABSTRACT 


18.  NUMBER  19a.  NAME  OF  RESPONSIBLE  PERSON 

Alejandro  A.  Hnilo 

PAGES  - - - 

19b.  TELEPHONE  NUMBER  (Include  area  code) 

005411  4709  8100  ext.1322 


DISTRIBUTION  A:  Distribution  approved  for  public  releat 


Standard  Form  298  (Rev.  8/98) 

Prescribed  by  ANSI  Std.  Z39.18 
Adobe  Professional  7.0 


INSTRUCTIONS  FOR  COMPLETING  SF  298 


1.  REPORT  DATE.  Full  publication  date,  including 
day,  month,  if  available.  Must  cite  at  least  the  year  and 
be  Year  2000  compliant,  e.g.  30-06-1998;  xx-06-1998; 
xx-xx-1998. 

2.  REPORT  TYPE.  State  the  type  of  report,  such  as 
final,  technical,  interim,  memorandum,  master's  thesis, 
progress,  quarterly,  research,  special,  group  study,  etc. 

3.  DATES  COVERED.  Indicate  the  time  during  which 
the  work  was  performed  and  the  report  was  written, 
e.g.,  Jun  1997  -  Jun  1998;  1-10  Jun  1996;  May  -  Nov 
1998;  Nov  1998. 

4.  TITLE.  Enter  title  and  subtitle  with  volume  number 
and  part  number,  if  applicable.  On  classified 
documents,  enter  the  title  classification  in  parentheses. 

5a.  CONTRACT  NUMBER.  Enter  all  contract  numbers 
as  they  appear  in  the  report,  e.g.  F33615-86-C-5169. 

5b.  GRANT  NUMBER.  Enter  all  grant  numbers  as 
they  appear  in  the  report,  e.g.  AFOSR-82-1234. 

5c.  PROGRAM  ELEMENT  NUMBER.  Enter  all 
program  element  numbers  as  they  appear  in  the  report, 
e.g.  61101A. 

5d.  PROJECT  NUMBER.  Enter  all  project  numbers  as 
they  appear  in  the  report,  e.g.  1F665702D1257;  ILIR. 

5e.  TASK  NUMBER.  Enter  all  task  numbers  as  they 
appear  in  the  report,  e.g.  05;  RF0330201;  T4112. 

5f.  WORK  UNIT  NUMBER.  Enter  all  work  unit 
numbers  as  they  appear  in  the  report,  e.g.  001; 
AFAPL30480105. 

6.  AUTHOR(S).  Enter  name(s)  of  person(s) 
responsible  for  writing  the  report,  performing  the 
research,  or  credited  with  the  content  of  the  report.  The 
form  of  entry  is  the  last  name,  first  name,  middle  initial, 
and  additional  qualifiers  separated  by  commas,  e.g. 
Smith,  Richard,  J,  Jr. 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND 
ADDRESS(ES).  Self-explanatory. 


8.  PERFORMING  ORGANIZATION  REPORT  NUMBER. 

Enter  all  unique  alphanumeric  report  numbers  assigned  by 
the  performing  organization,  e.g.  BRL-1234; 

AFW  L-TR-85-40 1 7-Vol-2 1  -PT-2 . 

9.  SPONSORING/MONITORING  AGENCY  NAME(S) 
AND  ADDRESS(ES).  Enter  the  name  and  address  of  the 
organization(s)  financially  responsible  for  and  monitoring 
the  work. 

10.  SPONSOR/MONITOR'S  ACRONYM(S).  Enter,  if 
available,  e.g.  BRL,  ARDEC,  NADC. 

11.  SPONSOR/MONITOR'S  REPORT  NUMBER(S). 

Enter  report  number  as  assigned  by  the  sponsoring/ 
monitoring  agency,  if  available,  e.g.  BRL-TR-829;  -215. 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT.  Use 

agency-mandated  availability  statements  to  indicate  the 
public  availability  or  distribution  limitations  of  the  report.  If 
additional  limitations/  restrictions  or  special  markings  are 
indicated,  follow  agency  authorization  procedures,  e.g. 
RD/FRD,  PROPIN,  ITAR,  etc.  Include  copyright 
information. 

13.  SUPPLEMENTARY  NOTES.  Enter  information  not 
included  elsewhere  such  as:  prepared  in  cooperation 
with;  translation  of;  report  supersedes;  old  edition  number, 
etc. 

14.  ABSTRACT.  A  brief  (approximately  200  words) 
factual  summary  of  the  most  significant  information. 

15.  SUBJECT  TERMS.  Key  words  or  phrases  identifying 
major  concepts  in  the  report. 

16.  SECURITY  CLASSIFICATION.  Enter  security 
classification  in  accordance  with  security  classification 
regulations,  e.g.  U,  C,  S,  etc.  If  this  form  contains 
classified  information,  stamp  classification  level  on  the  top 
and  bottom  of  this  page. 

17.  LIMITATION  OF  ABSTRACT.  This  block  must  be 
completed  to  assign  a  distribution  limitation  to  the  abstract. 
Enter  UU  (Unclassified  Unlimited)  or  SAR  (Same  as 
Report).  An  entry  in  this  block  is  necessary  if  the  abstract 
is  to  be  limited. 


DISTRIBUTION  A:  Distribution  approved  for  public  release 


Standard  Form  298  Back  (Rev.  8/98) 


FINAL  TECHNICAL  REPORT 

Contract/Grant  Title:  Nonlinear  dynamies  of  self-pulsing  all-solid-state  lasers. 
Contract/Grant  #:  FA9550-13-1-0120 
Reporting  Period:  March  2013  to  March  2015. 


Abstract. 

The  Project  involved  the  study,  both  experimental  and  theoretical,  of  extreme  events  (EE, 
sometimes  named  Optical  Rogue  Waves)  in  two  types  of  self-pulsing  lasers:  1)  all-solid-state 
(i.e.,  diode-pumped)  Nd:YV04+Cr:YAG  (“slow”  saturable  absorber);  and  2)  Kerr-lens-mode 
locked  Ti:Sapphire  (“fast”  saturable  absorber).  The  main  tool  is  the  analysis  of  the  time  series  of 
the  laser  pulses  and  of  the  images  of  the  laser  spots. 

Regarding  (1),  observations  on  a  specially  designed  and  constructed  prototype  showed  the 
existence  of  EE  in  the  pulse  energy  and  also  in  the  time  separation  between  pulses.  The  two 
types  of  EE  are  uncorrelated:  a  high  energy  pulse  is  not  preceded  by  a  longer  pumping  time. 
Therefore,  the  problem  is  to  identify  the  “reservoir”  for  the  energy-EE.  Correlations  were 
identified  between  the  spatial  complexity  of  the  laser  spot,  the  Eresnel  number  of  the  cavity,  the 
dimension  of  embedding  of  the  time  series  and  the  appearance  of  EE. 

Regarding  (2),  two  coexistent  modes  of  operation  are  observed:  transform-limited  pulses  (PI) 
and  chirped  pulses  (P2).  EE  are  observed  only  in  P2.  We  found  that  EE  arise  after  a  threshold 
similar  to  the  modulational  instability  is  crossed.  EE  exist  in  P2  only,  because  PI  is  unstable  for 
the  parameters’  values  above  that  instability  threshold.  If  the  system  is  forced  to  start  near  PI,  it 
evolves  into  P2  before  EE  can  be  observed  in  the  practice.  Einally,  the  observed  privileged 
distances  in  the  separation  between  successive  EE  are  the  residuals  of  the  “cold”  cavity 
periodicities,  perturbed  by  the  opposite  tendencies  of  an  expansive  Kerr  nonlinearity  and 
contractive  aperture  losses. 


Report  -  main  text  -  Introduction. 

The  subject  of  rogue  waves  (also  freak  waves)  has  recently  become  a  hot  topic  in  many  areas, 
including  such  dissimilar  fields  as  Optics,  Ecology  and  Stock  Exchange.  The  subject  was  bom  in 
Oceanography,  and  refers  to  the  observation  of  waves  of  large  amplitude  that,  although  rare, 
appear  more  often  than  one  would  expect  from  a  Gaussian  distribution.  “Optical  rogue  waves” 
were  defined  for  the  first  time,  as  large  fluctuations  in  the  light  intensity,  in  the  edge  of  the 
spectmm  produced  in  a  micro-stmctured  optical  fiber  pumped  with  femtosecond  laser  pulses  [1]. 

The  causes  and  mechanisms  of  formation  of  these  fluctuations,  in  Oceanography  as  in  the  other 
fields,  are  a  subject  of  research.  In  consequence,  it  is  unknown  whether  they  share  a  single 
description  or  not.  The  point  in  common  is  the  long  tailed  distribution.  In  Optics,  this  situation 
has  led  to  a  controversy  on  whether  the  term  “optical  rogue  wave”  is  appropriate,  or  not.  We  do 
not  want  to  take  a  position  in  this  controversy.  Hence,  we  will  use  the  neutral  term  “extreme 
event”  (EE)  in  this  Report.  Yet,  we  think  pertinent  to  note  that  the  Nonlinear  Schrodinger 
Equation  (NESE)  has  been  attempted  to  describe  the  formation  of  oceanic  rogue  waves  [2-3], 
and  that  the  same  equation  is  used  in  the  customary  approach  to  lasers  with  a  saturable  (i.e.. 
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nonlinear)  absorber  [4],  Therefore,  at  this  point  of  the  researeh,  it  should  not  be  disearded  the 
phenomena  to  have  in  eommon  more  than  it  seems,  espeeially  in  the  ease  of  the  Kerr-lens-mode- 
loeked  Ti:Sapphire  laser. 

In  this  Projeet,  we  studied  the  EE  observed  in  two  types  of  solid-state  lasers: 

1)  The  all-solid-state  (i.e.,  diode-pumped)  Nd:YV04  laser,  passive  Q-switehed  with  Cr:YAG  as 
the  saturable  absorber  (SA).  Typieally,  it  emits  pulses  of  10-100  ns  duration  at  a  repetition  rate  of 
5-20  KHz,  at  a  wavelength  of  1064  nm.  This  deviee  is  small,  effieient,  eeonomie  and  robust,  and 
is  used  in  many  applieations.  Knowing  the  meehanism  of  formation  of  EE  in  this  system  may 
lead  to  eontrol  their  emission.  This  would  allow  having  laser  pulses  of  high  intensity  at  times  of 
interest,  without  having  to  seale  up  the  whole  deviee.  It  would  be  an  important  advantage  in 
some  uses  of  interest  for  the  Air  Eoree,  for  example,  in  laser  rangefinders  or  target  illuminators 
aboard  small  unmanned  flying  vehieles.  This  laser  eorresponds  to  the  ease  of  a  “slow”  SA. 

2)  The  Kerr-lens-mode-loeked  (KEM)  Ti: Sapphire  laser.  Typieally,  it  emits  pulses  of  20-200  fs 
duration  at  a  repetition  rate  50-200  MHz,  at  a  wavelength  in  the  800  nm  range.  It  is  the  most 
widespread  souree  of  ultrashort  laser  pulses  nowadays.  It  is  also  elose,  in  mathematieal  terms,  to 
the  problem  of  the  oeeanie  rogue  waves  (through  the  tentative  eommon  deseription  with  the 
NESE).  Inereasing  the  knowledge  on  the  dynamies  of  this  deviee  is  of  praetieal  interest  in  itself 
Besides,  it  is  hoped  that  the  results  of  the  observations  on  this  laser  (whieh  are  mueh  faster  and 
easier  to  reeord  than  in  the  high  seas,  not  to  mention  the  ehange  of  the  eontrol  parameters)  may 
provide  some  elues  to  the  solution  of  the  problem  of  the  oeeanie  rogue  waves.  This  laser 
eorresponds  to  the  ease  of  a  “fast”  SA. 

In  the  following  lines,  we  briefly  report  the  results  of  our  studies  performed  on  the  EE 
observed  in  these  two  types  of  lasers.  Also,  in  the  all-solid-state  laser  with  loss  modulation.  This 
aetivity  was  not  planned  originally,  but  it  appeared  eonvenient  as  a  way  to  understand  the  more 
eomplex  ease  of  the  all-solid-state  laser  +  SA.  Publieations  and  other  seientifie  eontributions 
presented  during  the  Projeet  are  listed  at  the  end.  They  provide  further  details  on  the  studies 
briefly  reported  here. 

1)  EE  in  the  Nd:YV04+Cr:YAG  laser  (more  details  in  publications  #1,  #4  and  #5). 

We  designed  and  built  a  diode-pumped  Nd:YV04+Cr:YAG  laser  prototype  with  a  large 
Eresnel  number.  The  eavity  has  a  V-shape  and  the  Cr:YAG  erystal  (whieh  is  the  SA)  is  plaeed  in 
the  arm  where  the  mode  diameter  has  the  largest  variation  (see  Eigure  1).  This  allows  adjusting 
the  saturation  parameter  by  simply  displaeing  the  Cr:YAG  eloser  or  farther  from  the  output 
mirror,  with  a  mierometer  translation  stage.  Several  Cr:YAG  erystals  with  different 
transmissions  were  used.  The  aetive  medium  Nd:YV04,  instead,  is  always  the  same:  3x3x1  mm^, 
1%  Nd  and  eoatings  appropriate  for  longitudinal  pumping  with  a  2W  CW  laser  diode  at  808  nm. 
In  general,  this  prototype  operates  in  the  self-pulsing  (passive  Q-switehed)  regime. 

Depending  on  the  position  of  the  Cr:YAG  inside  the  eavity,  uniform  Q-switehing,  period 
doubling,  bifureations’  easeades  and  ehaotie  regimes  are  observed.  In  the  ehaotie  regimes,  there 
are  situations  with  and  without  EE.  The  elassieal  approaeh  to  deseribe  this  system  [5]  prediets 
the  bifureation  easeade  and  even  the  ehaotie  regime,  but  not  the  EE.  Our  earliest  observations 
indieated  that  EE  arise  easier  in  this  system  if  the  Eresnel  number  is  large,  suggesting  that  EE  are 
due  to  multimode  dynamies. 
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We  recorded  time  series  of  the  pulse  energies  (with  and  without  EE)  with  fast  photodiodes  and 
a  digital  oscilloscope.  After  processing  the  series  to  take  into  account  only  the  maxima  of  the 
pulses,  we  calculated  their  dimension  of  embedding  (dE)  and  the  values  of  the  Eyapunov 
exponents  (this  was  not  possible  for  all  the  experimentally  recorded  series,  probably  due  to 
noise).  Chaotic  and  hyperchaotic  series  were  observed.  EE  were  prone  to  appear  in  series  where 
dE  was  larger  than  the  value  of  4  of  the  classical  approach  in  [5]. 


Figure  1:  Scheme  of  the  laser  prototype.  LD:  pump  laser  diode,  GL:  GRIN  lens,  ND:  Nd:YV04  active 
medium.  Ml:  spherical  mirror  (R=250  mm),  M2:  output  mirror  (plane),  SA:  Cr:YAG  crystal,  T=90%,  L: 
lens,  FP:  fast  photodiode,  CCD:  camera,  PM:  power  meter,  a=20°,  Ll=130  mm,  L2=75  mm.  The  position 
X  of  the  SA  is  adjustable  with  a  micrometer  translation  stage,  to  obtain  different  dynamical  regimes. 


We  drew  return  maps  for  the  energy  and  time  (i.e.,  the  separation  between  successive  pulses) 
variables,  and  heterodyne  interferograms  with  a  modified  Mach-Zehnder  and  a  CCD  camera.  We 
also  recorded  the  self-correlation  and  the  mutual  information  of  the  time  series  obtained  at 
different  points  in  the  output  beam.  The  purpose  was  to  determine  the  importance  of  the 
multimode  dynamics  in  the  formation  of  the  EE. 

We  observed  that  not  only  the  pulse  energy,  but  that  also  the  variable  “time  between  successive 
pulses”  can  be  chaotic,  and  that  it  can  also  display  EE  (i.e.,  there  are  extraordinarily  long  times 
between  pulses).  Noteworthy,  the  energy-EE  are  not  anticipated  by  a  time -EE.  This  implies  that 
there  is  some  unknown  mechanism  of  energy  storage.  As  energy-EE  appear  in  large  Eresnel 
number  configurations,  we  speculate  that  this  mechanism  is  related  with  regions  of  the  active 
medium  that  are  not  depopulated  by  the  average  pulses  due  to  the  destructive  interference  among 
transversal  modes  (building  in  this  way  a  Lippman  grating).  These  regions  would  be  effectively 
depopulated  by  the  EE  pulses.  This  idea  is  supported  by  the  observation  that  energy-EE  are 
followed  by  a  long  time  before  the  next  (non-EE)  pulse. 

The  regime  with  EE  shows  partial  transversal  decoherence  (Eig.2,  Eeft).  Besides,  the  self¬ 
correlation  of  the  time  series  obtained  at  different  points  in  the  laser  spot  decays  in  the  regimes 
with  EE.  This  indicates  that  EE  in  this  system  are  related  with  spatio-temporal  chaos.  Yet,  we 
also  observed  that  the  time  series  is  more  regular  and  predictable  in  the  neighborhood  of  an  EE 
than  an  average  pulse,  even  inside  the  same  time  series  (Eig.2,  Center  and  Right).  This  suggests 
that  the  process  of  formation  of  EE  follows  a  deterministic  dynamics,  and  that  there  is,  in 
consequence,  some  hope  to  predict  and,  eventually,  to  control  the  EE  in  this  system.  As  it  is 
mentioned  in  the  Introduction,  this  possibility  is  of  practical  interest. 
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Figure  2:  Left:  heterodyne  interferograms  of  the  output  beam  (the  laser  spot  without  interferenee  is  on  the 
left)  for  different  dynamical  regimes.  Up:  intermittent  chaos  with  no  EE,  the  interference  figure  covers  the 
whole  spot,  the  field  is  transversally  coherent.  Down:  chaos  with  EE,  the  region  surrounded  with  the 
black  line  does  not  display  interference  lines,  hence,  it  is  not  coherent  with  the  reference  beam.  Center: 
superposition  of  112  oscilloscope  traces  centered  in  each  of  the  EE  of  a  time  series  with  a  total  of  lO"* 
pulses.  The  horizontal  line  indicates  the  position  of  the  EE  threshold  for  this  series.  Right:  superposition 
of  1 12  oscilloscope  traces  centered  in  pulses  of  average  intensity,  extracted  from  the  same  time  series  as 
in  the  figure  at  the  center.  The  dynamics  becomes  unpredictable  much  faster  for  the  average  intensity 
pulses  than  for  the  EE. 

Due  to  the  time  resolution  of  our  eamera,  the  interferograms  as  the  ones  in  the  Fig.l  average 
the  images  formed  by  many  pulses.  Yet,  we  manage  to  know  from  the  energy-return  maps  that 
the  EE  are  not  assoeiated  with  a  determined,  fixed  spatial  pattern.  It  remains  to  know  if  these 
“EE -patterns”  are  uniformly  brighter  than  the  average  or,  instead,  if  they  inelude  “hot  spots”,  that 
inerease  the  total  reeorded  energy  (like,  say,  a  supernova  inereases  the  luminosity  of  a  whole 
galaxy).  We  foresee  obtaining  single-pulse  images  of  the  output  laser  beam  to  determine  which 
alternative  is  the  correct  one. 

This  part  of  the  Project  is  the  core  of  the  plan  for  the  Ph.D  in  Physics  (in  progress)  of  Carlos 
Bonazzola,  to  be  presented  at  the  Universidad  de  Buenos  Aires. 

2)  EE  in  the  TiiSapphire  laser  (more  details  in  publications  #2  and  #3). 

The  EE  in  this  laser  take  the  form  of  pulses  of  high  energy.  The  time  separation  between  pulses 
is  fixed,  and  given  by  the  laser  cavity’s  round  trip  time.  Our  measurements  were  performed  on  an 
owner-made  KEM  TiiSapphire  laser  of  standard  X-configuration  (see  Eigure  3),  CW  pumped  at 
532  nm  by  a  5W  diode-pumped  and  frequency  doubled  Nd;YV04  commercial  laser  (Verdi).  The 
main  control  parameter  is  the  group-velocity-dispersion  (GVD),  which  is  adjusted  by  displacing 
one  or  both  of  the  intracavity  prisms  perpendicular  to  the  cavity’s  axis.  This  laser  is  not  placed  in 
our  lab,  but  in  the  Laboratorio  de  Electrdnica  Cudntica  of  the  Universidad  de  Buenos  Aires. 

The  KEM  TiiSapphire  laser  is  known  to  display  three  coexistent  modes  of  operation; 
continuous  wave  output  (named  PO),  transform-limited  fs  pulses  (PI)  and  positive-chirped  pulses 
(P2).  It  is  possible  to  go  from  one  mode  to  the  other  by  slightly  perturbing  the  laser  (say,  by 
tapping  the  positioning  device  of  any  one  of  the  mirrors).  The  laser  also  evolves  spontaneously 
from  one  mode  to  the  other,  in  a  time  scale  of  minutes,  most  probably  due  to  mechanical  noise 
and  thermal  drifts.  The  modes  PI  and  P2  have  their  own  characteristic  road  to  chaos:  PI  through 
quasi-periodicity,  P2  through  intermittency.  Be  aware  that  “mode”  here  means  a  dynamical  way 
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to  operate,  and  not  a  spatial  oscillation  eigenstate,  as  the  usual  Gauss-Laguerre  modes.  The  laser 
is  controlled  here  to  oscillate  always  in  a  single  transversal,  nearly  TEMoo  mode. 


Figure  3:  Scheme  of  the  Ti:Sapphire  laser.  LB:  pump  focusing  lens  (^lOcm),  M2, 3  spherical  mirrors 
(R=10  cm)  MPi,2  plane  FIR  mirrors.  Pi, 2  pair  of  prisms.  Distances  (in  mm):  M3-R=R-M2=50,  M2- 
MP2=140,  MP2-Mi=  465,  M3-Pi=297,  Pi-MPi=198,  MPi-P2=415,  P2-M4=109.  The  prisms’  positions  are 
adjusted  to  get  negative  total  GVD. 

Our  group  pioneered  the  observation  of  EE  in  this  laser  (in  fact,  in  any  laser  of  standard 
design  operated  under  normal  conditions)  [6].  Two  intriguing  features  were  noted:  one,  the  EE 
are  observed  in  the  chaotic  regime  of  the  mode  P2  only.  Two,  the  time  separation  between 
successive  EE,  as  measured  in  numbers  of  intermediate  non-EE  pulses,  is  always  a  simple 
combination  of  the  numbers  1 1  and  12  (which  we  called  “magic  numbers”). 


•  Mode  Pj 

•  ModeP, 


Figure  4:  Left:  SPM  and  energy  of  each  pulse  in  a  chaotic  time  series  with  10^  mode-locking  pulses. 
Almost  all  the  pulses  in  the  mode  P2  (in  red)  have  a  larger  value  of  SPM  than  even  the  largest  in  P 1  (in 
black).  The  spiraling  orbit  corresponds  to  a  quasi-periodic  excursion  of  the  mode  PL  Right:  fast 
evolution  of  the  system  from  PI  to  P2  when  the  laser  parameters  are  above  the  “modulational- instability¬ 
like”  threshold.  The  transition  to  P2  occurs  after  218  round  trips  (»2.5  gs  in  real  time),  total  length  of  the 
run:  lO'^  pulses. 


During  this  Project,  by  using  a  well-tested  description  of  this  laser  with  a  five-variables 
iterative  map  [7],  we  were  able  to  explain  both  intriguing  features.  In  the  Figure  4  (Eeft)  the  self- 
phase-modulation  (SPM)  and  energy  of  each  pulse  in  a  series  of  10^  are  plotted,  for  each  mode  of 
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operation  (PI:  black,  P2:  red).  It  is  noticeable  the  spiraling  orbit  corresponding  to  a  quasi¬ 
periodical  excursion,  typical  of  the  PI  mode.  The  EE  threshold  for  the  mode  P2  is  plotted  as  a 
vertical  dotted  line.  There  are  331  red  points  (they  seem  fewer  because  of  the  scale  of  the  figure) 
to  the  right  of  this  line,  indicating  that  there  are  331  EE  in  this  series.  The  mode  PI,  as  said, 
displays  no  EE  and  is  much  less  spread  in  the  energy  variable.  Even  though  the  P2-pulses  are 
spread  both  in  SPM  and  energy,  the  EE  all  have  nearly  the  same  value  of  SPM  which,  in 
addition,  is  relatively  low  for  the  mode  P2  but  higher  than  any  pulse  in  the  mode  PI.  This  value 
of  SPM  extends  to  the  left,  inside  the  region  of  the  non-EE  pulses,  and  forms  a  sort  of  “limit 
line”.  We  hypothesize  that  this  value  of  SPM,  about  10'^  fs'^,  corresponds  to  some  kind  of 
threshold. 

In  general  terms,  the  Modulational  (or  Benjamin-Eeir)  instability  of  a  periodical  solution  of  the 
NESE  occurs  when  the  nonlinearity  overwhelms  the  dispersion  scaled  with  the  frequency  of  the 
periodical  stable  solution.  In  our  case,  the  nonlinearity  is  the  SPM  and  the  GVD  is  a  constant 
parameter  (the  stable  mode-locking  regime  is  the  periodical  solution),  so  that  it  is  reasonable  to 
expect  a  dynamical  threshold  associated  to  the  minimum  value  of  the  SPM  in  a  given  mode.  This 
analogy  is  important  for  the  purposes  of  this  Project,  because  the  Modulational  Instability  is 
suspected  to  be  involved  in  the  formation  of  oceanic  rogue  waves. 

Even  though  no  pulse  in  the  PI  mode  has  a  SPM  above  10'  fs'  in  the  Pig. 4,  some  of  them  are 
just  at  the  border.  In  principle  at  least,  a  small  increase  in  the  pulse  energy  may  allow  the  PI 
mode  to  cross  the  threshold.  We  therefore  increase  the  small  signal  gain  in  the  numerical 
simulation  for  the  PI -map  (to  increase  the  pulse  energy)  and  looked  for  EE.  In  fact,  a  40% 
increase  in  the  small  signal  gain  does  produce  EE  in  the  mode  PI  too.  Yet,  EE  are  not  observed 
in  the  mode  PI  in  the  practice,  because  PI  is  unstable  for  this  value  of  the  gain.  This  is  best  seen 
in  the  Pig.4  on  the  right,  where  the  numerically  obtained  trajectory  of  the  system  in  the  phase 
space  is  plotted  (the  variables  are:  U:  energy,  S:  inverse  of  the  spot  area,  x:  pulse  duration,  they 
are  all  scaled  to  their  values  at  the  fixed  point  of  PI).  Even  though  the  system  is  forced  to  start  in 
the  neighborhood  of  the  fixed  point  PI,  it  rapidly  evolves  into  P2.  Translated  to  the  real  time 
scale,  the  system  does  not  remain  in  the  mode  PI  more  than  3  ps,  too  short  to  be  observed.  By 
the  way,  note  the  vertical  excursions  (the  vertical  axis  corresponds  to  the  energy)  in  the  mode  P2, 
which  corresponds  to  EE,  and  that  the  crossing  from  one  mode  to  the  other  occurs  precisely  after 
a  high  energy  excursion  (the  early  stage  of  an  EE?)  in  PI. 

We  conclude  that  the  formation  of  EE  in  this  laser  is  related  with  the  crossing  of  a  threshold  in 
the  SPM  value,  which  has  similarities  with  the  Modulation  Instability  of  the  NESE. 

Regarding  the  “magic  numbers”,  we  develope  a  description  with  a  simplified  iterative  map, 
which  has  the  form  of  a  (perturbed)  Moebius  map  in  the  complex  plane.  The  analysis  of  the 
solutions  of  this  map,  and  its  comparison  with  the  real  laser’s  design  parameters,  show  that  the 
magic  numbers  correspond  to  the  geometrical  periodicities  [8]  of  the  empty  cavity  (or  “cold 
cavity”,  i.e.,  with  no  nonlinear  effect  taken  into  account).  Therefore,  the  observed  dynamics  are 
the  residuals  of  the  periodic  orbits  when  subjected  to  two  opposite  perturbations:  one  is 
dissipative  (or  contractive  of  the  phase  space),  due  to  the  unavoidable  presence  of  transversal 
apertures  inside  the  cavity.  The  other  perturbation  is  expansive,  due  to  the  self-focusing  in  an 
appropriate  cavity’s  design,  which  is  found  to  be  precisely  that  of  the  real  one.  Einally,  we  check 
that  if  the  cavity’s  design  is  changed,  the  values  of  the  magic  numbers  change  accordingly. 
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3)  EE  in  the  Nd.'YVOd  laser  with  loss  modulation  (more  details  in  draft  #6  at  the  end  of  this 
Report). 

The  dynamics  of  the  formation  of  EE  in  the  transverse  multimode  Nd;YV04+Cr;YAG  laser,  as 
it  was  explained  in  the  Section  (1),  are  complex.  We  then  proposed  to  study  the  formation  of  EE 
in  a  simpler  system,  in  order  to  get  some  clues  to  understand  the  complete  process.  The 
underlying  hypothesis  is  that  EE  are  due  to  the  occurrence  of  an  external  crisis.  In  few  words,  an 
external  crisis  is,  in  a  multistable  system,  the  collision  of  an  unstable  orbit  belonging  to  one 
solution  with  the  chaotic  attractor  belonging  to  another  one.  The  collision  produces  an  abrupt 
increase  of  the  volume  of  the  phase  space  that  can  be  reached  from  the  chaotic  solution.  Recent 
theoretical  developments  show  that  EE  should  appear  in  this  case,  for  parameters’  values 
attainable  with  a  solid  state  laser  with  modulated  losses  [9]. 

We  then  designed  and  constructed  a  laser  prototype  similar  to  the  one  sketched  in  the  Fig.l. 

-5 

The  active  medium  is  again  a  Nd;YV04  3x3x1  mm  ,  1%  doped  crystal  pumped  by  a  2W  CW 
laser  diode  at  808  nm,  and  a  V-shaped  cavity.  There  are  important  differences,  however,  with  the 
prototype  used  in  the  Section  (1)  (which,  by  the  way,  has  not  been  dismantled);  the  Fresnel 
number  is  adjusted  to  be  low  (the  beam  is  nearly  single  transversal  mode  ~1.1),  there  is  no 
SA  crystal,  and  an  electro-optical  modulator  (EOM)  is  inserted  into  the  cavity.  The  EOM 
introduces  a  polarization  rotation  proportional  to  the  voltage  Wmod  provided  by  a  signal  generator. 
The  Nd;YV04  active  medium  ensures  that  the  field  is  linearly  polarized,  so  that  the  rotation 
introduced  by  the  EOM  is  equivalent  to  a  change  in  the  cavity  losses.  The  EOM  is  driven  at  a 
sinusoidal  frequency  /  close  to  the  frequency  of  the  relaxation  oscillations,  which  is  measured 
here  to  be  111  ±  IKHz.  As  before,  the  laser  output  is  detected  with  a  fast  photodiode  and  time 
series  recorded  in  a  digital  oscilloscope.  The  amplitude  of  the  modulation  Vmod  (i-C-,  the  output 
voltage  from  the  signal  generator)  is  the  control  parameter. 


Figure  5:  Left:  A  section  of  the  oscilloscope  trace  recorded  for  Wmotn  1-7  V,  y=100  KFlz,  Pump  power= 
1.024  W.  The  triangles  indicate  the  position  of  the  maxima  of  the  pulses,  which  compose  the  time  series 
that  is  numerically  studied.  Right:  Number  of  EE  and  average  intensity  as  a  function  of  the  loss 
modulation  amplitude  near  the  crisis,/=100  KFlz,  P=1.024  W. 

Depending  on  the  value  of  V„,oa  the  laser  output  is  observed  to  be  CW,  periodically  pulsed, 
and  chaotically  pulsed  with  and  without  EE.  At  low  values  of  V,„o*  pulsed  regimes  of  different 
periodicities  are  observed  to  coexist.  Uniform  pulsing,  2-period,  3-period,  until  6-period  regimes 
are  stable.  Higher  8-  and  10-period  are  also  observed,  but  they  are  rather  unstable.  It  is  possible 
to  pass  from  one  pulsed  regime  to  the  other  by  blocking  and  unblocking  the  cavity,  or  by  tapping 
one  of  the  mirrors.  Transitions  among  the  different  regimes  occur  spontaneously  too.  These 
observations  confirm  that  the  system  is  multistable. 
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As  Vmorf  is  increased,  bifurcations  are  observed  and  the  output  beeomes  aperiodic.  As  before, 
the  peaks  of  eaeh  pulse  eompose  a  time  series  (see  Fig. 5,  Left)  whieh  is  numerieally  analyzed. 
For  a  typieal  time  series,  we  obtain  dE=5  (the  surrogate  series  has  no  definite  value  of  dE)  and 
the  Lyapunov  exponents:  0.32,  0.11,  -0.08,  -0.30,  -0.76.  The  regime  is  henee  dissipative  and 
(hyper)ehaotie,  but  no  EE  are  observed  (Wmod  =  L6V).  As  Nmod  is  inereased,  EE  begin  to  appear. 
The  number  of  EE  is  maximum  at  2.12  V,  then  deereases,  and  EE  disappear  at  2.3V,  being  the 
regime  always  ehaotic  (see  Eig.5,  Right).  This  is  preeisely  the  behavior  predieted  for  this  system 
as  it  passes  through  an  external  erisis  [9]. 

The  study  of  the  behavior  of  the  pulses  in  the  neighborhood  of  an  EE  (similar  to  the  one  on  the 
right  side  of  the  Eig.2  obtained  for  the  laser  +  SA)  shows  that,  also  in  this  ease,  EE  are  more 
regular  and  predietable  than  the  average  pulse.  This  suggests  that  the  evolution  of  EE  oeeurs  into 
a  well  defined  manifold  in  phase  space,  supporting  the  hypothesis  of  their  deterministic  nature, 
and  giving  some  hope  to  the  goal  of  predieting  and  eontrolling  them. 

In  eonelusion  for  this  Seetion,  the  laser  with  modulated  losses  appears  as  a  promising  “toy 
system”  to  get  insight  into  the  more  general  properties  of  the  proeess  of  formation  of  EE  in  all¬ 
solid-state  lasers. 
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We  present  experimental  data  supporting  the  hypothesis  that  extreme  events  (or  optical  rogue  waves)  in  a 
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larger  than  theoretically  expected.  We  compare  such  time  with  theoretical  and  experimental 
measurements  of  the  positive  Lyapunov  exponents.  Furthermore,  we  explore  theoretically  a  possible 
protocol  to  maximize  the  predictability  of  an  extreme  event. 
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During  the  last  years  several  optical  systems  have  been  used  as  test  benches  for  the  study  of  extreme 
and  rare  events  [1-4].  Such  events  have  been  also  called  optical  rogue  waves  for  their  high  amplitude, 
much  higher  than  expected  in  “normal  circumstances”.  Furthermore,  it  is  already  accepted  that  optical 
rogue  waves  may  appear  as  a  consequence  of  very  different  types  of  dynamics.  They  can  be  expected 
relatively  frequent  in  linear  systems,  in  some  sense  induced  by  random  boundary  conditions  and  the 
possibility  of  multiple  interferences  [5-6].  They  have  been  predicted  and  observed  in  nonlinear 
dynamical  systems,  for  example,  after  the  development  of  a  modulational  instability  [7-10].  In  other 
words,  extreme  events  (EE)  appear  in  complex  systems,  and  most  of  them  are  reasonably  studied  only 
through  statistical  methods  and  represented  by  an  infinite  number  of  variables.  However,  some  systems 
can  be  modeled  by  ordinary  differential  equations.  The  simplest  system  of  this  category  may  be  a  laser 
with  a  modulated  parameter  [11-14].  The  appearance  of  rogue  waves  in  such  system  was  studied 
theoretically  [14]  and  the  appearance  of  EE  was  associated  to  an  external  crisis  of  a  strange  attractor.  EE 
have  been  theoretically  predicted  to  exist  in  a  CO2  or  solid-state  laser  with  a  modulated  parameter  [13]. 
The  numerical  simulations  indicate  the  appearance  of  EE  to  be  caused  by  an  external  crisis,  i.e.,  the 
collision  of  an  unstable  orbit  belonging  to  one  of  the  solutions  of  the  multistable  system,  with  the  chaotic 
attractor  belonging  to  another  solution.  The  collision  produces  an  abrupt  increase  of  the  region  that  can  be 
reached  in  the  phase  space,  and  hence,  to  EE.  As  the  control  parameter  gets  farther  from  the  point  of 
collision  the  number  of  EE  decreases,  even  if  the  average  amplitude  of  the  pulses  increases.  In  this  paper, 
we  present  experimental  evidence  in  support  of  this  explanation  of  the  origin  of  EE  in  an  all-solid  state 
with  modulation  of  losses. 
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The  quantitative  definition  of  an  EE  usually  is:  a)  amplitude  higher  than  twiee  the  “significant  wave 
height”  or  “significant  intensity”  I;/j  [7,8,14],  which  is  the  average  calculated  among  the  set  of  the  1/3 
highest  events  in  the  series  [15].  An  event  is  then  considered  “extreme”  if  its  abnormality  index  AI  = 
leven/Iy/j  IS  larger  than  2.  Alternatively:  b)  amplitude  higher  than  four  times  the  standard  deviation.  These 
two  definitions  can  be  coincident  or  not,  depending  on  the  form  of  the  distribution.  In  this  paper,  we  use 
the  second  alternative  unless  specified  otherwise. 


Figure  1 :  Scheme  of  the  experiment. 

The  scheme  of  our  modulated  all-solid-state  laser  is  shown  in  the  Fig.l.  A  Nd:YV04  crystal  is 
pumped  with  power  P  with  a  CW  diode  laser  at  808  nm.  The  Nd:YV04  material  ensures  that  the  field 
oscillating  into  the  optical  cavity  is  linearly  polarized.  A  standard  V-shaped  cavity  produces  a  waist  near 
the  output  (plane)  mirror,  where  an  electro-optic  modulator  (EOM)  is  placed.  The  EOM  introduces  a 
polarization  rotation  proportional  to  the  voltage  Wmod  provided  by  the  signal  generator.  As  the  field  is 
linearly  polarized,  this  rotation  is  equivalent  to  a  change  in  the  cavity  losses.  The  EOM  is  driven  at  a 
sinusoidal  frequency  / near  the  frequency  of  the  relaxation  oscillations,  which  is  measured  here  to  be  111 
±  IKHz.  The  laser  output  is  detected  with  a  fast  photodiode  and  time  series  recorded  in  a  storage 
oscilloscope.  The  typical  output  power  is  200  mW  for  P=1.5W.  The  beam  has  a  uniform  transversal 
shape,  M^=  1.1  ±  0.05.  Depending  on  the  values  of  the  control  parameters  P,  Nmod  and/  the  laser  output 
is  observed  to  be  CW,  periodically  pulsed,  and  chaotically  pulsed  with  and  without  EE.  With  f=100  KHz 
and  low  values  of  Wmod,  pulsed  regimes  of  different  periodicities  are  observed  to  coexist.  Uniform  pulsing, 
2-period,  3-period,  until  6-period  regimes  are  stable.  Higher  8-  and  10-period  are  also  observed,  but  they 
are  rather  unstable.  It  is  possible  to  pass  from  one  pulsed  regime  to  the  other  by  blocking  and  unblocking 
the  cavity,  or  by  tapping  one  of  the  mirrors.  Transitions  among  the  different  regimes  occur  spontaneously 
too,  in  a  time  range  of  the  order  of  the  minute.  The  pulses  have  duration  between  0.8  and  1.2  ps.  These 
observations  mean  that  the  system  is  multistable.  In  this  respect,  this  system  behaves  the  same  as  the  CO2 
laser  with  modulated  losses  [16]. 

As  Wmod  is  increased,  bifurcations  are  observed  and  the  output  becomes  aperiodic.  The  zoom  of  a 
typical  aperiodic  oscilloscope  trace  is  displayed  in  the  Figure  2.  The  peaks  of  each  pulse  compose  a  time 
series.  For  the  complete  time  series  of  this  example  (which  is  made  of  89602  peaks),  the  value  of  the 
dimension  of  embedding  and  the  Eyapunov  exponents  are  calculated  to  be  dE=5  (the  surrogate  series  has 
no  definite  value  of  dE),  A,,=  0.32,  0.11,  -0.08,  -0.30,  -0.76.  The  regime  is  hence  dissipative  and 
hyperchaotic. 
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Figure  2:  A  section  of  the  chaotic  time  series  recorded  for  N mocr  1-7  V,  7=100  KFlz,  P=1.024  W.  The 
triangles  indicate  the  position  of  the  maxima  of  the  pulses,  which  compose  a  time  series. 


Inside  the  chaotic  regime,  EE  are  observed  for  some  values  of  the  parameters.  In  the  Figure  3, 
histograms  of  the  distribution  of  the  pulses’  amplitudes  are  shown  as  y mod  is  varied.  The  regime  is  always 
chaotic,  but  it  displays  no  EE  at  1.8V,  it  does  at  2.12  V  and  again  it  does  not  at  2.3V.  This  is  the  behavior 
predicted  for  this  system  near  an  external  crisis  [13]. 
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Figure  3:  Flistograms  of  the  pulses’  amplitudes  for  increasing  loss  modulation  amplitude.  The  dotted  line 
indicates  the  threshold  for  EE. 


In  the  Figure  4,  we  plot  the  number  of  observed  EE  as  y mod  is  scanned  through  the  crisis.  Compare 
with  the  Figure  8  in  the  Reference  13.  As  predicted,  the  number  of  EE  increases  as  the  crisis  is 
approached,  and  decays  after  the  system  has  climbed  to  the  next  “step”  of  higher  intensity.  Note  that  the 
average  output  intensity  increases  monotonically. 

In  a  numerical  simulation  of  the  semiconductor  laser  with  optical  feedback  in  the  short-cavity  regime 
[17],  EE  are  calculated  to  be  more  predictable  than  the  average  chaotic  pulse.  The  same  feature  is 
observed  here.  In  the  upper  part  of  the  Figure  5,  the  traces  corresponding  to  the  74  EE  of  the  time  series 
corresponding  to  the  peak  in  the  Fig.4  are  displayed.  The  pulses  neighboring  the  EE  occur  at  nearly  the 
same  time  and  have  nearly  the  same  intensity.  Naturally,  they  “blur”  as  the  distance  to  the  EE  increases. 
For  the  average  pulse  instead  (lower  part  of  Fig.5)  the  characteristics  of  the  neighboring  pulses 
immediately  become  unpredictable. 
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Figure  4:  Number  of  EE  and  average  intensity  as  a  function  of  loss  modulation  amplitude  near  the  crisis, 
7^100  KHz,  P= 1.024  W. 
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Figure  5:  Superposition  of  74  EE  and  the  same  number  of  pulses  of  intensity  near  the  average  intensity, 
both  extraeted  from  the  same  time  series,  which  corresponds  to  the  maximum  in  the  Fig.4. 

This  suggests  that  the  evolution  of  EE  occurs  into  a  well  defined  manifold  in  phase  space.  This  result 
supports  the  hypothesis  of  a  deterministic  nature  of  the  observed  EE,  and  give  some  hope  to  the  goal  of 
predicting  and  controlling  them. 

The  study  of  the  EE  in  the  laser  with  a  modulated  parameter,  as  a  subject  of  dynamical  systems,  is 
interesting  by  itself  Besides,  it  will  lead  to  a  deeper  understanding  of  the  operation  of  these  lasers  and, 
eventually,  to  an  improvement  of  their  performance.  This  simple  system  may  also  provide  a  convenient 
test  bench  to  study  the  general  features  of  the  phenomenon  of  EE.  For,  lasers  dynamics  evolve 
instantaneously  (at  the  human’s  timescale)  and  their  control  parameters  are  easy  to  adjust. 
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appearance  of  EE. 

Regarding  (2),  two  coexistent  modes  of  operation  are  observed:  transform-limited  pulses  (PI)  and  chirped 

pulses  (P2).  EE  are  observed  only  in  P2.  We  found  that  EE  arise  after  a  threshold  similar  to  the 

modulational  instability  is  crossed.  EE  exist  in  P2  only,  because  PI  is  unstable  for  the  parameters’  values 

above  that  instability  threshold.  If  the  system  is  forced  to  start  near  PI ,  it  evolves  into  P2  before  EE  can  be 
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EE  are  the  residuals  of  the  “cold”  cavity  periodicities,  perturbed  by  the  opposite  tendencies  of  an  expansive 
Kerr  nonlinearity  and  contractive  aperture  losses. 
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